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Abstract 

We analyze dissipative scale effects within a one-dimensional theory, devel¬ 
oped in [L. Anand et al. (2005) J. Mech. Phys. Solids 53], which describes 
plastic flow in a thin strip undergoing simple shear. We give a variational 
characterization of the yield (shear) stress — the threshold for the inset of 
plastic flow — and we use this characterization, together with results from [M. 

Amar et al. (2011) J. Math. Anal. Appl. 397], to obtain an explicit relation 
between the yield stress and the height of the strip. The relation we obtain 
confirms that thinner specimens are stronger, in the sense that they display 
higher yield stress. 
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1 Introduction 


A number of experiments have shown that conventional plasticity fails to capture 
the size-dependent behavior of metallic specimens undergoing plastic flow in the size 
range below 100 microns, with smaller samples being, in general, stronger (see [20] 
for a review). 

Substantial theoretical work has been carried out to extend conventional plasticity 
to the micron scale. It is acknowledged that size effects observed in metallic samples 
are associated to the inhomogeneity of plastic flow [5], a fact that motivates a number 
of strain-gradient plasticity theories, starting with Ref. [lO] . 

In the so-called non-local or high-order theories, the flow rule that governs the 
evolution of plastic strain is a partial differential equation which requires the speci- 
hcation of appropriate boundary conditions. The hrst of such theories was proposed 
by Aifantis [1]; the vast majority of subsequent high-order theories were derived using 
the virtual-power principle, by taking into account power expenditure by higher-order 
stresses that are work-conjugate to the plastic-strain gradient [SI [H EZl ESI DS]. 

Apparently, the theories developed by Gurtin and Anand [ISl US] are those that 
have inspired most mathematical work. One of its distinctive aspects is that the 
full plastic distortion (the sum of a symmetric plastic strain and a skew-symmetric 
plastic spin) is accounted for. In particular, the issues of existence and uniqueness 
of solutions for strain-gradient plasticity with plastic spin, as considered in Ref. [I8] . 
has been addressed in Ref. [8] in the case of two-dimensional setting of anti-plane 
shear, and in Refs. mi.izTi, and [2B] in the full three-dimensional setting. The model 
for plastically-irrotational materials proposed in Ref. [I9| was studied in Ref. m 
Of particular importance for the present paper are the existence theorems for strain- 
gradient plasticity based on the notion of energetic solution, which have been proved 
both in the small-strain [T6| and in the large-strain [23] setting. 

The flow rules proposed in Ref. [T8l |19| are of particular interest because they 
incorporate two length scales: 

• an energetic scale L, which appears from letting the free-energy density depend 
on derivatives of the plastic strain, E^, through the Burgers tensor, G = curlE^; 

• a dissipative scale i, which arises from letting the gradient of plastic strain 
rate, VEp, enter the dissipation-rate density. 

The form of the free energy density is motivated by dislocation mechanics. In 
particular, the choice of letting the free energy to depend on plastic strain gra¬ 
dient through the Burgers tensor follows from the presumption that the so-called 
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geometrically-necessary dislocations (whose density is measnred by G) play a major 
role in determining size-dependent response, a presnmption that hnds its jnstihcation 
in homogenization resnlts from discrete-dislocation models [I51E3]. 

Becanse of the complicated natnre of the non-local flow rnle, it is not easy to 
nnderstand how its solntions are affected by the material scales. On the other hand, 
snch nnderstanding is crncial in order to identify these scales by comparison with 
experiments. Thus, parallel with the literature dealing with modeling, researchers 
have also endeavored to investigate how the various scales may affect the nature of 
solutions, not only for the Gurtin-Anand theory, but also for other strain-gradient 
plasticity theories. 

This task is usually accomplished by working out a simple analytical problem 
that mimics some experimental setup. For example, scale dependence for the torsion 
experiment was investigated in Ref. [22] (by numerical and asymptotic considera¬ 
tions) in the framework of the Fleck & Willis theory [T3| and in Ref. [9] (by rigorous 
arguments) for energetic scale effects within the Gurtin-Anand theory [T9|. More¬ 
over, for the distortion-gradient plasticity (which accounts also for plastic spin), 
specihc hnite-element schemes for the torsion problem have been recently proposed 
in Ref. [7| . Problems involving micro-bending have been scrutinized in Ref. |2I] and, 
more recently, in Ref. [13] in the case of non-monotone loading. 

With a similar goal in mind, a simplihed flow rule, formulated in one spatial 
dimension, was derived and analyzed in Ref. [3] to investigate the effects of both 
the energetic and the dissipative scales. Such flow rule, which mimics the traction 
problem in simple shear symmetry, will be introduced in Section In the same sec¬ 
tion we also make a comparison with conventional plasticity. From our comparison 
two facts emerge: 1) that the length-scale £ is expected to be a source of additional 
strengthening; 2) that the natural way to quantify strengthening is to consider in¬ 
crease of the Yield stress, Ty, i.e., the value of the (shear) stress that triggers plastic 
flow in an initially virgin sample. 

The aim of this paper is to rigorously prove that, according to the flow rule we 
consider, Ty is strictly increasing with i, that is to say, smaller samples are stronger. 
In fact, within this simplihed framework of simple shear symmetry, we will be able 
to determine the dependence of Ty on i explicitly. Results and proof are stated (in 
renormalized variables) in Section ^ which also contains an outline of the arguments 
(details are given in Sections EjlS ). In summary, we will hrst argue that Ty may 
characterized as the smallest value that the renormalized dissipation 
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attains among all 0 G +h)) snch that 0(?/)d?/ = 1 (see Theorem 1). This 

constrained minimization problem had already been introdnced in |1] and analyzed 
in [2], showing that a minimnm is attained in BV, which is smooth in the interior 
and satishes the corresponding Enler-Lagrange eqnation. We will then argne that 
these resnlts permit to explicitly characterize ry in terms of I (see Theorem and 
Fignrej^. 


2 Problem setup 

2.1 The traction problem 

The one dimensional theory developed in Ref. |1] describes plastic flow in a body 
having the shape of an inhnite strip of width 2h, namely, = |x = {x, y, z) G : 
-h <y < h], as sketched in Fig. 1. We restrict attention to the so-called traction 



Fignre 1: An inhnite strip of height 2h, clamped on the bottom side and snbject to 
a uniform shear traction r on the top side. 

problem^ describing an ideal experiment in which the bottom side of the strip is 
clamped and a uniform shear stress r is prescribed on the upper side. We work in 
the rate-independent setting of quasistatic evolution in plasticity and we limit our 
attention to the case of proportional loading, that is to say, we assume that r be 
strictly increasing with respect to time. With this assumption, we may label each 
instant by the corresponding value of the shear stress and adopt r in place of time 
as the independent variable. 

Because of translational invariance in the x- and ^-directions, it is natural to look 
for solutions that enjoy the same invariance properties. Precisely, we assume that 
the two kinematic helds of interest, namely the displacement u and the plastie strain 
EP, have the following representation: 

u = u{y, r)ei, EP = '-f^{y, r)sym(ei 0 es). 
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with {ej : i = 1,2, 3} the canonical basis of Consistent with this assnmption, we 
take the stress tensor T to be spatially constant, and having the representation 

T(r) = r(ei ® e 2 + 62 ® ei). 


2.2 A local flow rule: strengthening and hardening 

If the material is modeled in the framework of Mises plasticity with kinematic hard¬ 
ening, the flow rnle governing the evolntion of the shear strain 7 ^ may be written 
as 

f r - 

\^GSign(7P), 

where So > 0 is the coarse-grain yield strength, n is the kinematic-hardening coeffi¬ 
cient, a snperimposed dot denotes differentiation with respect to the loading param¬ 
eter r, and 

r {+ 1 } ifx> 0 , 

Sign(x) = < [—1, -1-1] if a: = 0, 

[ {— 1 } if a: < 0 . 

Granted that the body is in its virgin state at the beginning of the experiment, 
namely, 

7P(2/,0) = 0, (2) 

the solntion of Q is easily worked ont and, on introdncing the positive-part operator 
(■)+= max{-, 0 }, can be written as 




{t/Sq - 1)^ 

K 


This solntion displays the typical featnres of a stress-strain diagram from classical 
plasticity; in particnlar: 


• the increase of S'o is associated to strengthening, that is, an increase of the 
threshold for the inset of plastic flow, the Field shear stress: 

ty = Sq, (3) 

• the increase of n, with Sq hxed, is associated to hardening, that is, an increase 
of the shear stress reqnired to attain a given amonnt of plastic shear. 


5 



( 4 ) 


On multiplying ([^ by 7 P, we obtain the free energy balance 

is„K('|;(77) +Sol7'’l=Tr, 

the free energy density being given by The balance (|^ can thus be in¬ 

terpreted as a splitting of the internal power ryP expended on plastic flow into an 
energetic part and a dissipative part^ T-dis^yp _ 5*0 lyP]. Accordingly, we may say that, 
in the present context. 


• strengthening is a dissipative effect, whereas 

• hardening is an energetic effect. 


2.3 A non-local flow rule: size-dependent strengthening and 
hardening 

Using the strain-gradient plasticity theory of Ref. [Ij we derive in Appendix 1 a 
non-local, rate-independent flow rule. In particular, we replace the hrst of ([^ with; 


T 


- So {Kjf - = t"”' - kf‘, 

and the inclusion in ([^ with: 


Sn 


e Vers (yP, £yP) , 


(5a) 


(5b) 


where the index y denotes partial differentiation with respect to y and 


Vers(v) = /{H} if|v|#0. 

[ {v G ; |v| < 1} if |v| = 0. 

Problem (|^ must be complemented by both initial conditions, for which we again 
choose the virgin-state condition ([^, 

Y\t=o = 0, (6a) 


and boundary conditions, for which we choose microscopic hard conditions: 


'y^\y=—h 'y'^\y=-\-h 0 * 


(6b) 
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As explained in Appendix 1 , the partial differential equation (5a) is constitutively 


augmented microforce balance engendered by a version of the principle of virtual 
powers that accounts for power expenditure on the time derivative of the shear-strain 
gradient 7 ^. In particular, taking the formal variation of the plastic free energy 


E'-en = f 




— h 


+ L‘^ (7y)) dy 


(7) 


and dehning the plastic dissipation rate 




>i»‘(7'‘) = So / +em?)^dy, 


( 8 ) 


' —h 


the following identity is arrived at: 




^p(^p) + ^P(^P) = 
dr 


T^dy, 


( 9 ) 


' —h 


which is again interpreted as a splitting of work expenditure (the right-hand side of 
([^) into an energetic part and a dissipative part. Given that L, resp. i, appear 
in the energetic, resp. dissipative part of the energy balance ([^, in line with the 
discussion in 92.2 


one may expect that the extra energy brought into play by L enhances harden¬ 
ing effects, and that the extra dissipation associated to £ is a source of additional 
strengthening. 


The role of L has been recently scrutinized in [9], rigorously conhrming this 
expectation in the case of torsion of thin wires. The role of i has been investigated 
both formally and numerically in |1]. In view of the discussion in (2.2 (cf. in 
particular a natural way to rigorously quantify the role of i is to determine how 
the yield shear stress. 


ty ■= sup {r > 0 : 7 ‘’(-, r) = O}, (10) 

depends on £. This constitute the goal of this paper. 
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Main results 


3.1 Scaling 

In order to single ont the relevant parameters, we introdnce dimensionless indepen¬ 
dent variables: 

= ^ ft = — 
h’ 5, • 


0 


Consistent with this choice, we introduce the dimensionless parameters: 




( 11 ) 


The nonlocal flow rule ([^ can now be reformulated in the domain / := (—1, -f-l) and 
takes the form (henceforth, for typographical convenience, we drop the superscript 
p from yP): 

+ A\r = , 

(7,7r). 


( 12 ) 


where the index r denotes partial differentiation with respect to r. Initial and mi¬ 
croscopically hard boundary conditions ([^ now read as 


7(r,0) =7(-l,«) =7(+l,«) = 0 (r,») 6 / x |0,+oo) 


( 13 ) 


and the renormalized plastic free energy, resp. dissipation-rate, are given by 

B(7) := I / (7^ + Ay?) dr, I”(7) := / 7/7^ + V7jdr ( 14 ) 

(cf. 0, resp. (|^). In renormalized varibales, our aim becomes that of rigor¬ 
ously quantifying the dependence on the renormalized dissipative scale. A, of the 
renormalized yield shear stress (cf. @) 

^ — By sup {(1 > 0 : 7(1!) — 0 G [ 0 , 11 ]}, ( 15 ) 

namely, the largest value attained by the renormalized shear stress 6 prior to the inset 
of plastic flow. 



3.2 Energetic formulation 

We assume hereafter that k > 0, A > 0, and A > 0. Being a rate-independent 


dynamical system, the flow rule (12)-(13) can be formulated in many equivalent 


ways. The formulation that best suits our needs is the so-called energetic formulation 
proposed in Ref. |2S]- With a view towards formulating (12)-(13) in the energetic 


format, we introduce the (renormalized) energy functional: 

S’{9,'y) := E{')) — 0 [ ydr. 


(16) 


As usual, we write 7 ( 0 ) := 7 ( 0 ,-). We can now give the definition of energetic 
solution. 

Definition 1 (Energetic solution). Given 0 > 0, a function 7 : [0,0] —)■ Hq{I) is 


an energetic solution to (12)-(13) if the function [0,0] 39 h->■ ^( 6 *, 7 ( 6 *)) = — Jjjdr 


is in T^((0, 0)) and if the following conditions are satisfied for all 9 G [0, 0].- 

^{9,'y{9))<^{9,v) + df{-f{9)-v) for all v e (17a) 

^(0,7(0)) + disif(7;[O,0]) = - r [j{'9)drd^, (17b) 


where disii,( 7 ; [ 0 , 0 ]) is the total variation of 'y on [ 0 , 0 ] with respect to the distance 
6^(71,72) = ^(71 - 72)1 he., 


N 


disvi>( 7 ; [ 0 , 0]) := sup ^(7(0^) - 7(0j-i)) : A^ G N, 0 = 0o < • • ■ < 07v = 0 

In the present setting (quadratic energy) the next proposition is established with¬ 
out burden by invoking, for instance. Theorem 2.1 in Ref. 123- 


Proposition 1. There exists a unigue solution 7 of (12)-(13). Moreover, 9 h->■ 7 ( 0 ) 
is Lipschitz continuous as a function from [0, 0] to 

3.3 Characterizations of ry 

The first main result of this paper is the following characterization of 0y: 


(15). Then 


Theorem 1. Let 7 be the unigue energetic solution to (12)-(13) and let 9 y as in 
9y = inf |t(0): 0 G jfidr = l| . (18) 
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In order to explain the relation between the two quantities, it is convenient to 
briefly illustrate the main steps in the proof, whose details are given in ® We begin 


by observing that the energy-balance condition (17b) is identically satished for all 


Q e (O,0y). Thus, what determines the inset of plastic flow is the loss of stability of 
the trivial state 7 = 0 . This leads us to consider the stability indicator-. 


m{9) -.= inf where <h 0 ((/)) := ^( 6 *, 0 ) + \ 1 '( 0 ). 


(19) 


We will indeed argue that 


6y = inf {0 > 0 : m{6) < 0 } 

(cf. Proposition]^. Next, we note that the plastic dissipation rate T is (positively) 
homogeneous of degree one in 7 , whereas the plastic free energy E is quadratic. Then, 
a simple scaling argument can be used to show that the reduced stability indicator 


m{9)-.= inf where <h 6 )((/)) ;= T((/)) — 0 (f)dr 

Ji 


( 20 ) 


is equivalent to the stability indicator-. 

m(r) < 0 77 fh(r) < 0 

(cf. Proposition]^. The last step of our argument consists in observing that, again 
by homogeneity, for negative values of m we can restrict our attention to the subspace 

of tests (j) satisfying the normalization condition J (pdr = 1 : this leads to Theorem 

/ 


3.4 The formula for Ty 

The second main result of this paper is the following explicit 


Theorem 2. The renormalized yield shear stress 9 y 


diassipative scale A 


£ 

h 


are related by 


IL 

So 


formula for ry: 
and the renormalized 


A 


71 { 9 y — 


^y5ZL 

■\/6*y — 1 ) + 29y arctan 



( 21 ) 
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The proof is provided in Section]^ and relies on results from Ref. [2], guaranteeing 
that the relaxation in BV (/) of the inhmum problem in (18) admits a minimizer (py 
which is smooth in I and satishes the Euler-Lagrange equation 


9y = 


4>y 


X 


dr 




( 22 ) 


By a suitable change of dependent variable, we converts ( |22[ ) into a first-order dif¬ 
ferential equation with two side conditions. The extra side condition selects the 


eigenvalue Oy of the E-L equation (22), yielding (22). 


The graph of ry/S'o, recovered from (21), is plotted in Fig. (recall (0 and 
([T^). Our result conhrms that as the sample becomes smaller, i.e. X = i/h increases, 
the actual yield strength increases: hence smaller samples are stronger. Needless to 
say, the results from out plot agree with the numerical calculations carried out in 
Ref. [1] and reported in Figure 4 thereof. 



Figure 2: Solid line: renormalized effective yield strength ty/Sq versus renormalized 
dissipative scale £/h, as from formula (21). Dashed line: the upper bound ^ < 1-|- | 
derived in [5]. This plot agrees with the result computed numerically in [5] and 
reported in Fig. 4 thereof. When comparing the two hgures, the reader should take 
into account that in the present paper the symbol h denotes half the thickness of the 
strip, whereas in [5] the same symbol denotes the overall thickness. 
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Our explicit formula provides additional insight concerning the asymptotic be¬ 
havior of the actual yield strength for small and large values of h. In particular, from 
(21) one hnds that, for 0 < — 1 -C 1, 

A 


~ -^\/Oy — 1 , 

TT 


which implies that, for 0 < A -C 1, the renormalized actual yield strength has the 
following asymptotic behavior: 


TT 


6 'y — 1 ~ 0 < A <C 1. 

We also note that, as X = i/h ^ oo, a. linear relation is recovered: 
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TT 


Y 


A ~ for A 3> 1. 
4 


4 Proof of Theorem [U 


Existence and uniqueness of the minimum in (19) is readily ascertained through the 


direct method of the calculus of variations, owing to coercivity, lower semicontinuity, 
and convexity of <h 6 i in Hq{I): 

Lemma 1. For any A > 0 there exists a unique minimizer (px of the infimum problem 


in (19). 


The hrst step is to show that if the trivial state is not stable at a certain value of 
the renormalized shear stress 6 during the loading process, then it is not stable for 
whatever higher value: 


Lemma 2. The function [0, 0] 36 i—)■ m{6) defined in (19) is non-increasing. 


Proof. Let (fe be the unique minimizer of <l>e. First, we observe that 

(J)/)> 0 a.e. in / for all 6 > 0 . 


(23) 


Indeed, obviously 0o = 0; for 6 ^ > 0, if 06i < 0 in a set J of positive measure, then 


(by the dehnitions (16) and (14) of 8 , resp. T) we would have <h6i(|06i|) < ^e{,4>e), a 


contradiction. Thus, given 9i < 62 , we have 


m{92) - 

^e2{4>e2) 


< 

*^ 6 l 2 ( 06 »i) 

(by def. of (pg^) 



f n \ 
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as desired. □ 

The previous lemma is expedient to arrive to the following characterization of By- 


Proposition 2. Let 7 be the unique energetic solution to (12)-(13) and let By and 


m as in (15), resp. (19). Then 

By = inf {0 > 0 : m{B) < 0} . 


Proof. Let us set 6 ^ = inf {0 > 0 : m{B) < 0}. We notice that, since m{B) is non¬ 
increasing, m{B) = 0 in [0,0). Hence, by direct substitution into 0 , we see that 
the trivial function 0 i—)- 0 is an energetic solution on the interval [0,0). By the 


uniqueness of the energetic solution, and by (15), it follows that By > 0. 


The reverse inequality follows from the monotonicity of 0 h-)■ m{B): suppose 
indeed that 0 < By] then, by Lemma there exists B < By such that m(0) < 0; 


however, B < By implies that 7 ( 0 ) = 0; thus, by (17a) and (19), this means that 
m( 0 ) = 0 , whence a contradiction. □ 


We now show that the reduced stability indicator dehned in (20) can be used to 


detect the inset of plastic flow. Indeed, we have the following equivalence: 


Proposition 3. 


By = inf{0 > 0 : m{B) < 0}. 


(24) 


Proof. In view of Proposition it suffices to show that 

m( 0 ) <0 if and only if m( 0 ) < 0 . 

Since by dehnition $0 < <I)e, m(0) < 0 obviously implies m{B) < 0. For the reverse 
implication, let us assume m(0) < 0. Then there exists (f G Hq{I) such that $ 0 ( 0 ) < 
0. On the other hand, by the 1-homogeneity of $ 0 , 

, $ 0 ( 0 ( 5 ) ~ 

hm -= < 0 . 

a-^o-h a 

Thus ^g{a(j)) < 0 for o; > 0 sufficiently small, whence m{B) <0. □ 

With Proposition at hand we are now ready to establish the variational char¬ 
acterization we have been after. 
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Proof of Theorem Let 


0y(A):=inf ^(0): (f e Udr = 1 


(25) 


On recalling the definitions of 'L and given in (14), respectively (20), we see that 
the ineqnality 

Oy < ^y(A) (26) 

is implied by the following chain of implications: 


9y{\) < 9 


4/(0) <9 for some 0 e Hq{I) snch that / 0dr = 1 


inf 


6'0dr + 4/(0) ) < 0 


^ m{9) < 0 

9y ^ 9. 


Having established (26), it remains for us to prove the reverse inequality: 

9y ^ ^y(A)- 


To this aim, let 9 G (O,0y(A)). By the definition (25) of 0y(A), we have 
9 / 0dr = 0 < T(0) for all 0 G Hq{I) such that / 0dr = 1 . 


(27) 


(28) 


Since both sides of the inequality in (28) are positively 1-homogeneous, (28) upgrades 
to 

9 I 0dr < 4/(0) for all 0 G Hq{I) such that / 0dr > 0. (29) 


In turn, since T is non-negative, (29) upgrades to 

0 < ^(0) - 0 ^0dr ® $e(0) for all 0 G H^il) 

which holds for all 9 G (0, (A)). Summing up, we have the implication: 

0 < 6 ^ < 9y{X) rh{9) = inf ^e{4>) > 0 0 < ^y. 


(30) 


whence (27), since 6*y > 0 by definition. 


□ 
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5 Proof of Theorem [2] 


The infimum problem in (18) was addressed in Ref. |2]. 


Consider the relaxation of 


T(0) := inf | liminf : {0*.} C lTo’^(/), 0 in T^(/)|, (31) 

l. fc—^oo ) 

i.e. the largest lower semicontinous extension of T. It is shown in Ref. [2] that the 
relaxation T has the following representation for 0 G i?R(/)|^ 

= I + dr + A||r>*^||(/) + A (|.^.(-1)| + |.^.(+1)|). (32) 


Notice that, as is customary in the BV setting, homogeneous boundary conditions 
are now incorporated in the functional through the penalization term X\(j)\{dl) = 
A|0 — 0\{dl), which measures the jump between the trace of 0 and the prescribed 
null value. 

The following results were established in Ref. [2]. 


Theorem 3 (see Thm. 5.1 in Ref. 12]). Let T as in ( [3l| . 
G SBV{I) such that J^0ydr = 1 and 


There exists a unique 


T(0y) = min |tI/(0) : 0 G L^{I), 

Moreover, (py is even, strictly decreasing in [0,1), and smooth in (—1,1); further¬ 
more, it solves the Euler-Lagrange equation 

I 1 

= -A^f- , * (33) 

and it satisfies 



lim Mb 

r^l- 0y(O) 


Oy - 1 


e 


Y 


and 


lim 



— (X). 


(34) 


^Here ||;A|j denotes the total variation of a measure g (see e.g. 0 Def. 1.4]) and resp. 
denote the absolutely continuous, resp. singular, part of Df with respect to the Lebesgue measure 
(see e.g. [31 Th. 1.28 and §3.9]). We also refer to [3] for definitions and basic properties of the 
spaces BV{I) and SBV{T). 
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Remark 1. Notably, (34) shows that the solution 0y G SBV{I) of the relaxed 
minimization problem does not satisfy the boundary conditions 0(—1) = 0(1) = 0; 
generally speaking, this amounts to saying that, in order to minimize ^ with mass 
constraint, paying a jump discontinuity at the boundary is cheaper than attaining 
the boundary value zero. 

We are now ready to prove Theorem 

Proof of Theorem^ In view of Theoremand since Hq{I) is dense in BV{I), 

T(0y) = 9y- (35) 

We also notice that, since d0y/dr < 0 in [0,1) and 0y is positive with fj 0y(r)dr = 1, 


Oy ® ^(0y) 5 ^ 


r = 1. 


(36) 


Now, consider the function 


d<^y 

dr 


C(r):=-A-p 

Since 0y is smooth and positive in I, ( is smooth as well. We note that 






2 • 


Hence, since 0y > 0, 


<pY 




= WT 


(37) 


By making also use of the Euler-Lagrange equation, we see that C, satishes the fol¬ 
lowing differential equation: 


dr 


(j)Y 


It follows from (36) and (38) that ^ > 0. Hence, 

rC{i-) 




( 38 ) 


1 

A “ 


dc 

dr 


dr = 


- \/l - C J C{o) 9y - a/1 - 


dC- 
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In addition, since 0y is even and because of (34), we have that 


C( 0 ) = 0 and lim ({r) = 1 . 
1 —^ 1 — 


Therefore, 


1 

A 


dC 




( 39 ) 


The integral on the right-hand side of (39) is well defined and can be computed 
explicitly. As a result, we arrive at formula ( 21 ) for the renormalized actual yield 
stress. □ 


A The nonlocal flow rule 

In this section we briefly recapitulate the steps leading to the flow rule ([^, as devised 
in Ref. |1], with a few changes from the original path. At variance with the previous 
sections, we do not assume proportional loading. Accordingly, the independent vari¬ 
ables are now y G (—h, +h) and t, which stands for time, and the index y denotes 
partial differentiation with respect to y. 

A.l Principle of virtual powers 

We start from the decomposition 

Uy = 'f'' + yP (40) 

of the shear strain Uy into an elastic part 7 ® and a plastic part 7 P. This decomposition 
is accompanied by the prescription that, given any part P = {a,b) C {—h,+h), the 
internal power expended within P has the form: 

^nt(^) = J ry® -F rPyP -F k^^^dy. (41) 

Thus, power expenditure by the macroscopic shear stress r is accompanied by work¬ 
ing of the plastic microstress and gradient microstress kP. If body forces are left 
out of the picture, the external power expended on P = (a, h) is localized on the 
boundary dP = {a, b} and has the form: 

^ext(iP) = T{b)u{b) + k^{b)Y{b) - T{a)u{a) -k'^{a)Y{a), 

where r and fcP are, respectively, the macroscopic and the microscopic shear tractions. 
The application of the principle of virtual powers yields: 
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1 ) the identification between stress and traction 

T = T, 

along with the macroscopic-force balance: 

T, = 0; (42) 

2 ) the identihcation of kP with fcP, along with the microscopic force-balance: 

T = — kfy. 


A.2 Constitutive prescriptions 


Consistent with the choice (41) for the internal power expenditure, it is assumed in 
Ref. in that the free-energy density (p depends on the triplet ( 7 ®, 7 ^, 7 p through a 
constitutive equation of the form: 


^ = ^(7^7^,7^)• 

It is also assumed that the constitutive mapping delivering the free-energy density 
is the sum: 

^(7^7^7P) = ^(7'^) + ^"(7^7P) 

of an elastic-energy mapping which takes into account the elastic shear, and a 
defect-energy mapping fP, which depends on the plastic shear and on its gradient. In 
particular, the elastic-energy mapping is given the form ^( 7 ®) = ^^( 7 ®)^, with G > 
0 the shear modulus. This assumption is accompanied by the standard constitutive 
prescription r = whence: 

r = G7^ (43) 

The microstresses are then split into an energetic part and a dissipative part by 
setting 

rP = + r“, fcP = + k^'^, 

where 

en _ r en _ 

^ ayP ’ 

so that the following reduced form of the dissipation inequality is arrived at: 

0 < 
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By analogy with the constitutive equations describing viscoplastic behavior in metals, 
the following constitutive equations have been considered in Ref. |1]: 


dP 






5 = H{S)dP, 


5(0) = 5o. 


dP\^il 


(44) 


Here: S is the current yield strength, an internal variable whose value at time f = 0 
is equal to the initial yield strength So and whose time derivative is proportional to 
the effective flow rate dP through a (isotropic) hardening/softening function H{S); 
do is the reference flow rate; m > 0 is the rate-sensitivity parameter. 

The constitutive prescription (5b) follows by setting H{S) = 0 (no isotropic 
hardening) and by formally letting m —)• 0 in (44) (rate-independent limit). The 
partial differential equation (|^) is recovered by choosing 


A.3 The traction problem 

In the traction problem, the bottom side of the strip is clamped, that is, 

u{—h, t) = 0, (45) 

and a time-dependent shear traction Thif) is prescribed on the upper side, that is, 

T{h,f) = Th{t). 


On recalling that the shear stress is spatially constant by (42), we see that the shear 
stress r(f) appearing in the flow rule ([^ is a prescribed, spatially-constant held. 
Thus, the how rule (44) can be solved for the plastic shear yP without knowing the 
displacement held. The latter is recovered by integrating (40) and (43), and by 
taking (45) into account, that is to say. 


u{y,t) = 


' —h 


G 


+ 7 ^(s,f) ) ds. 


Acknowledgements. G.T. thanks Lallit Anand, Lorenzo Bardella, and Patrizio 
Neh for stimulating discussions on strain-gradient plasticity. 
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